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We investigate the doping of a geometrically frustrated spin ladder with static holes by a com- 
plementary approach using exact diagonalization and quantum dimers. Results for thermodynamic 
properties, the singlet density of states, the hole-binding energy and the spin correlations will be 
presented. For the undoped systems the ground state is non-degenerate, with translationally invari- 
ant nearest-neighbor spin correlations. In the doped case, we find that static holes polarize their 
vicinity by a localization of singlets in order to reduce the frustration. This polarization induces 
short range repulsive forces between two holes and an oscillatory behavior of the long range two- hole 
energy. For most quantities investigated, we find very good agreement between the quantum dimer 
approach and the results from exact diagonalization. 

PACS numbers: 75.10.Jm, 75.30.Hx, 75.40.Mg 
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I. INTRODUCTION 

Geometric frustration is a key factor, leading to exotic 
phases in quantum spin systems. Valence bond (VB) or- 
dering, with and without breaking of discrete lattice sym- 
metries, is one possible type of ground state symmetry. 
Examples are the spin-1/2 zig-zag ladder—, the checker- 
board lattice^ii^, the ji-j2-j3 model^i^, or the Shastry- 
Sutherland model^. Another exotic phase is the spin Hq- 
uid (SL) which displays no apparent magnetic order but 
may break topological symmetries, which presumably is 
the case for the two-dimensional (2D) spin-1/2 kagome 
antiferromagnet (AFM}ii^i^. 

Static holes, i.e. nonmagnetic defects, are an im- 
portant probe of such quantum phases. In case of en- 
hanced VB correlations or VB order without broken lat- 
tice symmetries, the situation is similar to spin-ladders 
and the 2D Heisenberg AFM, where static holes gener- 
ate spin-1/2 moments in their vicinit y^°i^^i^^ . The lat- 
ter induce triplet excitations at low-energies. In a SL, 
e.g. in the kagome AFM, a very different scenario has 
been observed, in which singlet pairs accumulate close 
to the holes, rather than spin-1/2 momentsi^. In case 
of VB states with spontaneous breaking of the lattice 
symmetrjii^, and other valence-bond ground states,— 
bound spin-1/2 moments also seem to be absent. Apart 
from the properties of a single static hole, correla- 
tion effects between two holes may shed Hght on the 
(de) confinement of mobile carriers in frustrated quantum 
magnetsi^ii^iii^'i^ - at least for kinetic energies small com- 
pared to the exchange coupling, i.e., t <^ J . Along this 
line, hole-hole interactions have been studied recently in 
the kagome AFM^^ and in quantum dimer models^^. In 
both cases, deconfinement of the static holes has been 
found. For two recent /iSR-studies of doped kagome 
AFMs, see refs. [Uii. 

Additional insight into the response of frustrated spin 
models to static holes can be obtained from a reduced 




Figure 1: ID Bow-Tie lattice structure of corner sharing tri- 
angles. Spin-1/2 moments reside on all vertices. Thick grey 
(black) bonds correspond to one QD configuration. Combi- 
nation of grey and black QDs represent one transition graph. 
Circles at the center correspond to two static holes. 

description of the singlet sector in terms of short-range 
resonating valence-bonds^^, so-called quantum dimers 
(QDs}^'2i. Classical dimer models on bipartite lattices 
allow for a height representation and tend to be confin- 
ing, with a Hnearly or logarithmically attractive 'string 
potential' between hole a^^i^^ . Non-bipartite lattices lack 
a height representation and allow for both, confining as 
well as deconfining22i^i22ii^i^ phases. Similar behavior 
is Hkely for QD models. 

In this work we intend to shed further light on the 
role of static holes in geometrically frustrated magnets by 
studying the spin-1/2 AFM on a bow-tie ladder (BTL). 
The BTL is a one-dimensional array of corner sharing tri- 
angles, as depicted in Fig. [1] and may be viewed as the 
medial lattice of the two-leg laddei>2^. The AFM Heisen- 
berg model on the BTL with nearest-neighbor exchange 
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has been investigated extensively by linear spin wave the- 
ory, exact diagonalization (ED) and density matrix renor- 
malization group (DMRG) in ref. [H. This analysis has 
uncovered several issues which remain open. In particu- 
lar, ED on up to iV = 30 sites suggests that the BTL is 
qualitatively similar to the Kagome AFM, i.e. the system 
is a SL with a spin gap and a large number of singlets 
(oc N^) below the first triplet. However, ED does not 
allow extrapolation of the triplet gap to the thermody- 
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Figure 2: Ground state energy versus inverse system size for 
iV = 12 - 30 (ED) and iV = 12 - 36 (QD). 

namic limit. DMRG up to iV = 120 indicates a vanishing 
spin gap for — > oo. Therefore, apart from a SL state it 
remains possible that the BTL is critical or displays VB 
order with a very large unit cell. 

Here, we will not focus on these properties of the clean 
BTL in the thermodynamic limit. Instead, we will inves- 
tigate the local correlations of two static holes - holons 
hereafter - introduced into finite BTLs. To this end, re- 
sults from two complementary methods will be discussed, 
namely, complete ED on the one hand, and restricted 
diagonalization in the QD subspace on the other hand. 
Recently, QDs on a BTL lattice have been considered, for 
the so-called /i-model of shortest length resonance-moves 
for the QDs. There, for certain ad-hoc resonance ampli- 
tudes the QD model was shown to be identical to the 
transverse-field Ising chain at criticality22. In contrast 
to that study, we will account for all resonance moves 
resulting from eqn. ([1]). 



II. THERMODYNAMIC PROPERTIES 

First, we briefly consider the clean system and com- 
pare the ground state energy and specific heat, obtained 
by ED, with results from the QD approximation. The 
generalized eigenvalue problem in the QD Hilbert space 
uses the transition graphs depicted in Fig. [T]and is set up 
according to ref. |33| (to which we refer the reader for de- 
tails). Periodic boundary conditions (PBC) apply to all 
results discussed in the following. For all systems studied 
we find the undoped ground state to be non-degenerate, 
with translationally invariant nearest-neighbor spin cor- 
relations, which is consistent with a spin liquid state or 
a valence bond crystal with a unit cell large compared 
to the system's sizes. In Fig. [21 we show the ground 
state energy Cg versus the inverse system size 1/A^. For 
QDs, proper N needs to be a multiple of 6. From Fig. [21 
a deviation of approximately 4-5% between ED and QD 
can be extrapolated in the thermodynamic limit. This 
has to be contrasted against the relative dimension D 
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Figure 3: Specific heat Cv versus temperature T. (a) and 
(b): comparison of ED and QD for = 12 and 18. (c) QD 
for N = 24-36. 



of the complete Hilbert space, Dq of the singlet sector, 
and Dqd of the QD space, which are D = 2^ , Dq = 
Nl/[{N/2y.{l + N/2)l], and Dqd = 2^/^+1^ p^^. 7v=30, 
this implies D/Dq/Dqd « 1.110V(9.7 10'')/(2.0 10^). 
In view of these ratios, the results in Fig. [21 are in reason- 
able agreement. The figure shows a very small 'even/odd' 
oscillation with respect to iV/6 in the QD results. 

A fingerprint of SL states in strongly frustrated mag- 
nets is the accumulation of singlet states at low energies. 
In the kagome AFM, the numer of singlets below the first 
triplet has been observed to grow exponentially with sys- 
tem size^. In the BTL, power-law behavior with an expo- 
nent 2 has been suggested^. The singlet accumulation 
at low-energies leads to a characteristic singlet-triplet 
double-peak structure in the specific heat Cy^ i^^i^'^i^^ . 
For the kagome AFM a substantial fraction of the singlet- 
peak in Cv is believed to be due to QD fiuctuationa^. 
For the BTL, this is an open issue. Therefore, in Fig. [3] 
we consider Cy from QD and ED calculations. The lat- 
ter agree with identical results from ref. [H. In view of 
large finite-size effects our emphasis is on a comparison at 
identical system sizes, rather than in the thermodynamic 
limit. First, there is a remarkable similarity between the 
positions of peaks in ED and QD in Fig. [3Ua,b). Second, 
for the leftmost peak in Fig. [3jb), at low temperatures, 
the magnitude of Cy and the entropy from ED and QD 
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Figure 4: Number of singlets in ED versus energy at = 18, 
contrasting holons (solid line) with 2 holons (dashed). Solid 
(dashed) vertical line: triplet gap for (2) holons. Holons are 
located at two adjacent sites on one BTL leg. Eo refers to 
ground state energy. 

are comparable. I.e., we conclude that the low-T peak in 
Cy of the BTL is caused primarily by QD fluctuations. 
As can be seen from Fig. [DJc), the low-T peak stabi- 
lizes as N increases. Third, in the high-energy region 
T > 0.5 a single peak can be observed which contains 
most of the entropy (note the log T-scale) . This entropy 
is significantly larger in the ED as compared to the QD. 
This difference is due to triplet excitations not present 
in the QD space. Fig. [3jc) shows that Cy from QD 
has converged to the thermodynamic limit for T > 0.1 
at -/V > 36 - a system size we cannot reach with ED. 
Finally, the intermediate-T feature in ED seems strongly 
system-size dependent, as already noted in ref. [H. 



III. STATIC HOLES 

In this section, we consider the change in the singlet 
density-of-states (DOS), the binding energy, and the spin 
correlations which stem from introducing two holons into 
the BTL. 




Figure 5: Number of QD singlets versus energy at = 18, 
contrasting holons (solid line) with 2 holons (dashed). 
Holons are located at sites identical to Fig. ID Eo refers 
to ground state energy. 

low-energy singlet DOS decreases, while the triplet gap is 
set by excitations distant from the holons and therefore 
less affected. While similar effects have been reported for 
the kagome AFMi^, this behavior is drastically different 
from VB states on bipartite lattices, where static holes 
tend to bind spin-1/2 moments which leads to an increase 
of low-energy triplet DOS^. 

While the QD variational space does not contain triplet 
excitations, it is nevertheless instructive to compare the 
impact of holons on the singlet DOS between the ED and 
QD calculations. In Fig. \5\ we show the number of sin- 
glets above the ground state energy in the QD case for 
a system size and a holon placement identical to that of 
Fig. m First, the overall trend for QDs depicted is re- 
markably similar to that in ED. Second, the suppression 
sets in at somewhat larger energies, of E-Eq ~ 0.3 for 
QDs, as compared to E-Eq ~ 0.2 for ED. Finally, it can 
be seen that in the energy window depicted the QD sin- 
glets make up for approximately 50% of the total number 
of singlets in ED. 



B. Two- holon Energy 



A. Density of States 

Fig. |4] shows the integrated singlet DOS versus energy, 
contrasting the case of zero holons against that of two 
holons placed at adjacent sites on the leg of a BTL with 
an even number of sites. The main point of this figure is 
the number of singlets below the first triplet. Evidently 
this number is reduced by introducing the holons. A 
similar effect is observed for any other relative separation 
of the holons and also in the case of a single holon. This 
low-energy suppression of singlet DOS originates from 
the reduction of frustration at the neighboring sites of 
the holons, which favors binding of singlets in the vicinity 
of the holons (see also section BlI Cp . Consequently the 



Now, we turn to the binding energies 

Ae^, (L) = ES (L) - ES {L = 1) (2) 

of two holons separated by L sites along the center (leg) 
= a of the BTL, where Eq refers to the corresponding 
ground state energies. In Figs. [6] and [7| we present two 
aspects of Ae2/j {L). In the former, we compare the bind- 
ing energies obtained from ED with those from QDs for 
various systems sizes. First, and remarkably, the ED and 
QD results are very similar, both for a = center and leg. 
Second, these figures demonstrate that the holons expe- 
rience a short range repulsive force, with a maximum 
binding energy at approximately 3 lattice sites. This is 
in sharp contrast to quantum AFMs with VB ground 
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Figure 6: Comparing two-holon energies at holon separation 
L from ED (solid line with open squares) with QDs (solid line 
with solid squares) for different systems sizes A'^=18, 24, and 
30. (a)-(c) [(d)-(f)] refer to two holons in [on] the center [leg] 
of the BTL. 

states on bipartite lattices. There, two holons experience 
a maximum binding energy if placed on those nearest- 
neighbor sites which are occupied by singlet dimers in 
the undoped case. Moreover, separating the holons will 
generate a string of 'misplaced' dimers of length propor- 
tional to the holon separation. This generates a confin- 
ing potential. Such confinement cannot be inferred from 
Fig. |6l Rather, from the largest system which allows for 
a direct comparison between ED and QDs, i.e. N = 30, 
an oscillatory behavior at large distances can be antici- 
pated. This can be corroborated by extending the QD 
calculation to larger system sizes as shown in Fig. |7l 
There, clear oscillations in Acg^ (L) can be observed in 
the right panel, both on the legs and in the center. In the 
next section, an intuitive picture of the two-holon ener- 
gies in terms of the dimer coverings will arise. From the 
preceding, it is tempting to speculate that mobile holons 
on the BTL are bound only weakly and will deconfined 
already for kinetic energies small compared to J. 

Finally, Fig. |7| points towards some of the limitations 
of the QD approach. While the qualitative behavior is 
identical for ED and QDs, i.e. panels (a) vs. (c) and 
(b) vs. (d), panel (d) shows a significant variation of the 
absolute value of the QD binding energies, depending on 
whether N/6 is even or odd. The latter even/odd effect 
is absent in the ED results. 



C. Spin-Spin Correlations 

The nearest-neighbor spin-spin correlations sim = 
(Si • Sm) are a direct measure of the singlet amplitude 
on the bond Im. Here, we consider the impact of two 
holons on s;™ at zero temperature. Fig. [8] summarizes 
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Figure 7: Evolution of two-holon energies at holon separation 
L with system size for ED (solid line with open squares) with 
A*'=18-30 in (a), (b) and for QDs (solid line with solid squares) 
with iV=18-42 in (c), (d). 

our results, both for ED and QDs on the largest system 
for which we have performed ED, i.e. N=30. In this fig- 
ure Sim is visualized in terms of 'bond-thickness'. First, 
we note that sim fulfills the sum rule 



Sir. 



(3) 



where Eq [L) refers to the ground state energy with two 
holons in a relative configuration denoted by L. Using 
the results from Figs. |6l|7l and|2l the sum rule has been 
checked to hold for all cases we have studied, both for ED 
and QDs and also for the undoped case, which will not be 
discussed here. From Fig. |8l it is evident that the holons 
introduce a strong polarization into the magnetic back- 
ground, which leads to spatial oscillations of the spin-spin 
correlations, with a pattern depending on the separation 
of the holons. This has to be contrasted against the un- 
doped case, in which sim is almost homogeneous along 
the BTL with only a very small transverse difference be- 
tween the central rungs and the legs. Specifically, the 
figure demonstrates that sim tends to be largest in the 
vicinity of the holons. This corroborates the picture of a 
doping induced reduction of singlet fluctuations, due to 
singlet-binding to the holons as discussed in Sec. IIII Al 
In fact, the maximal values of \sim\ in the vicinity of the 
holons as shown in the flgure are close to 3/4, implying 
a singlet amplitude of 1 on those bonds. Note that for 
a small number of bonds, ferromagnetic correlations of a 
rather small absolute magnitude arise, both in ED and 
for QDs. 

The tendency of the system to accommodate singlets 
next to the holons provides for a direct interpretation 
of the holon-repulsion, namely the number of 'strong' 
singlets which can be bound to the vicinity of the holons 
is less if their separation is small. This can be read off 
directly from Fig. [8] in going from panel 1) to 5). 
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a) b) 

Figure 8: Nearest neighbor spin-spin correlations sim from ED 
and QDs, for two holons (full circles) on the center and leg of a 
BTL with = 30 for a) and b). The width of solid lines (with 
vertical slash) is a linear measure of —sim{sim) for sim < 0(> 

0) . 1). . .5) labels all non-equivalent two-holon positions. For 
each 1). . .5) the adjacent lower (upper) graph corresponds to 
one ED (QD) result, both for a) and b). Numbers refer to 
maximum values of —sim observed for ED and QDs in a) and 
b). 

For the majority of holon placements depicted in Figs. 
IHK) and b), namely a) 1), 2), 3), and 5) as well as b) 

1) , 2), and 3), the qualitative real-space structure of sim 
is remarkably similar for ED and QDs. For group a) 4) 
there are some differences observable right of the second 
holon. In groups b) 4) and 5) there are clear qualitative 
differences between ED and QDs, especially on the leg 
free of holons. This suggests that QDs provide for a 
better description of holons in the 'bulk' of the BTL than 
on its boundaries. 

For a quantitative comparison between ED and QDs, 
we turn to Fig. [9l where we show several cuts along 
the BTL through the legs and rungs displaying the val- 
ues of —Sim versus a measure of distance given by a 'ID 
center-of-mass' CM of each bond. CM runs in steps of 
2 (1/2) from (-1/2) to 2A^/3 - 2 (2A^/3 - 3/2) on the 
legs (rungs) . The top four panels in Fig. [9] are represen- 
tative of those cases in Fig. [8] which suggest qualitative 
agreement between ED and QD. It is readily apparent 
that there is also excellent quantitative agreement. The 
lower two panels in Fig. [9] refer to —sim on the upper leg 
(rung) for the worst case of Fig. [8]b) 5), where ED and 
QD show strong qualitative differences. 

In the following, we provide additional discussion of 
the results shown in Figs. [6] - [H When both holons 
are positioned in the center axis (see Fig. [SK), reflec- 
tion symmetry around this axis is retained. This allows 
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Figure 9: Negative nearest neighbor spin-spin correlations 
—Sim along the upper(lower) leg(rung) for u(l)l(r) from 
ED(QDs) for open(closed) circles with two holons and A''=30, 
corresponding to Fig. [8]and panels a) 3), b) 1), and b) 5). CM 
labels 'center of mass' of bonds. Crosses mark bonds linked 
to holons where Sim has been set to 0. 

for the formation of highly symmetric and simple spin 
structures. First, as noted above, isolated singlets form 
at the bonds situated closest to the holons for all holon 
separations. Second, more subtle geometrical constraints 
decide if other spin structures (e.g., loops with even num- 
bers of spins or short chains) will form, where they will 
be located, and how much singlet correlations there will 
be between them. For example, in Fig. [8^) 3), at a holon 
distance of L = 3, apart from the four isolated dimers, 
one can clearly notice the presence of 4- and 8-spin loops 
in the shape of three diamonds and a four-pointed star, 
respectively. These structures are only weakly connected 
by singlet correlations^!. Obviously, the distance be- 
tween the holons is one of the controlling parameters in 
this respect. For example, at distance L = 5, open 3- 
spin chains are present in the shape of boomerangs, be- 
sides diamonds and dimers. However, in this case, these 
structures are not as independent from each other, i.e., 
the singlet links between them are stronger, as compared 
to L = 3. Therefore, some holon distances are more 
favorable than others in allowing for the formation of 
such independent 'frustration releasing' structures. This 
leads to the oscillating behavior of the binding energy 
in Fig. [T]:)^. Finally, since weakly connected dimers, 
boomerangs, and diamonds can be described rather well 
by the QD basis, the very good quantitative agreement 
between ED and QD for holons located in the center axis 
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is no surprise. 

When the holons are both on the same leg, reflection 
symmetry around the central axis is lost. Therefore, 
apart from the dimers close to the holons, simple struc- 
tures as described above are less favorable and form only 
for some holon distances, as for example the diamonds 
seen at L = 2 in Fig. [8b). Instead, and in contrast 
to the centered holons, more complex correlations of the 
spins occur, which involve the presence of 'spin chains' of 
considerable length. These cannot form when the holons 
are in the center axis because of symmetry considera- 
tions and avoidance of uncorrelated free spins. These 
chain structures with longer range spin-spin correlations 
are not well described by the QD basis. Their presence 
leads to quantitative discrepancies with the ED results. 
The difference in agreement between Fig. Elja) and (c) 
versus Fig. Eljb) and (d), as well as the discrepancies 
in the two lower two panels in Fig. [9] corroborate this 
interpretation^. 

IV. CONCLUSIONS 

In summary, we have performed a complementary nu- 
merical analysis of a geometrically frustrated quantum 
spin ladder with and without static holes using exact di- 
agonalization and a truncated basis of quantum dimers. 



For the undoped system, we have shown that dimers al- 
low for a reasonable approximation of the ground state 
energy and the low temperature speciflc heat. In the 
doped case, we have shown that the system can re- 
lease frustration through binding of singlets and other 
extended unfrustrated spin structures to the static holes. 
Results for the holon binding energies and the nearest- 
neighbor spin correlations confirm this picture. Analysis 
of dynamic correlation functions, as e.g. magnetic Ra- 
man scattering, and the finite size scaHng of overlaps be- 
tween ED and QD states, will be presented elsewhere^. 
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